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Abstract 

We evaluate the influence of thermal fluctuations on superconducting rings that enclose a 
magnetic flux, using the time-dependent Ginzburg-Landau (TDGL) or the Kramer-Watts- Tobin 
(KWT) model, while thermal fluctuations are accounted for by means of Langevin terms. This 
method is applicable in situations where previous methods are not, such as nonuniform loops, 
rings with large width to radius ratio and loops with large coherence length to perimeter ratio. We 
evaluate persistent currents, position and statistical behavior of flux-induced vortices and lifetime 
of metastable fluxoid states. The influence of nonuniformity on the persistent current does not 
depend strongly on the details of the cross-section profile; it depends mainly on its first harmonic, 
but not only on it. As a consequence of nonuniformity the maximum of the persistent current 
shifts to smaller fluxes and the passage between fluxoid states remains non-hysteretic down to 
lower temperatures than in the case of a uniform sample. Our results using TDGL agree remark- 
ably well with recent measurements of the persistent current in superconducting rings and with 
measurements of the position of a vortex that mediates between fluxoid states in an asymmetric 
disk with a hole; they could also provide a plausible explanation for unexpectedly short measured 
lifetimes of metastable states. Comparison between TDGL and KWT indicates that they lead to 
the same results for the persistent current, whereas KWT leads to larger lifetimes than TDGL. 

PACS numbers: 74.40. +k, 05.10.Gg, 73.23.Ra, 74.78.Na 
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I. INTRODUCTION 



Superconducting rings have attracted the interest of physicists during many decades, 
because they are an easily accessible system in which measurable quantities depend on 
the enclosed flux, rather than on local fields only. More precisely, these quantities depend 
on the magnetic flux modulo $0 = /ic/2e, as expected from quantum behavior. One of 
the milestones in this endeavor was the Little-Parks experiment, J. in which the transition 
temperature is an oscillatory function of the applied flux. 

A closer look at the Little-Parks results shows that, as expected from a one-dimensional 
system, there is actually no phase transition. The transition region is dominated by thermal 
fluctuations and what experimentalists call "transition temperature" is in most cases the 
temperature at which the resistance becomes some given fraction of the normal state resis- 
tance. Since the Little-Parks effect is most pronounced for small samples, fluctuations are 
especially important. Superconducting rings are a compelling system for studying fluctua- 
tions: theoretically, they are a nontrivial system which is not invariant under time reversal, 
has a rich yet simple phase diagram, and both Gaussian and non-Gaussian fluctuations may 
be involved; experimentally, the influence of fluctuations is quite directly controllable and 
detectable. 

A theory for the evaluation of fluctuation contribution to conductivity above the criti- 
cal temperature is due to Aslamazov and Larkin.-S An exact method for evaluation of the 
average current around a uniform loop, as a function of the temperature and enclosed flux, 
was developed by von Oppen and Riedel (vOR).- Although the vOR method is exact from 
the point of view of statistical mechanics (provided that the energy spectrum is described 
by the Ginzburg-Landau model), it is limited to a static situation and to rings that are 
perfectly uniform and one-dimensional; moreover, for small values of the parameter 7 (will 
be defined in Sec. Ill AD numeric implementation of this method becomes exponentially dif- 
ficult. An independent theory, which discusses the influence of the shape of the loop and 
uses a two-level approximation is due to Daumens et al.— Recently, the persistent current 
in a superconducting loop was analyzed by Schwiete and Oreg.- For the regime in which 
the superconductor-normal transition occurs at finite temperature, they find asymptotic ex- 
pressions for the current near integer or half-integer flux. For sufficiently small rings, there 
exists a quantum critical point at which the superconductor-normal transition occurs at 
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zero temperature; in this case they evaluate the persistent current near the critical point by 
means of the full fluctuation propagator. The current thus obtained is larger than predicted 
by classical fluctuations alone. 

The vOR theory was tested by Zhang and Price. ^ They measured the susceptibility and 
found values larger by an order of magnitude than predicted. Recently, Koshnick et ai— 
repeated this experiment for several samples with widths of the order of 100 nm, using a 
scanning SQUID, and found good agreement with the theory. Rosario et ai— measured 
conductivity in ultrathin cylinders at half-integer number of fiuxoids, and found exponential 
dependence on the temperature above the transition, in contrast with the power dependence 
predicted by Aslamazov-Larkin. 

In this article we study the influence of the width and nonuniformity of the loop on the 
average current and on the transitions between fluxoid states. We focus on the Little-Parks 
temperature range, i.e. the range of temperatures (depending on the flux) for which the 
superconducting transition would occur in the absence of fluctuations. For higher tem- 
peratures, currents are small and experimental and numeric values are noisy; for lower 
temperatures, there is hysteresis and the statistical average becomes irrelevant. Somewhat 
lower temperatures will also be considered, because fluctuations and nonuniformity lower 
the temperature at which hysteresis appears and, sporadically, we will enter the hysteretic 
range. Most of the article is devoted to evaluation of the current around the rings studied in 
Ref. and comparison with their experimental values. For most samples, our results agree 
with the vOR theory (and therefore with the theoretical analysis in Ref. Q); when they do 
not, we attribute the difference to the width and/or nonuniformity of the samples. 

The most widespread theoretical tool in the study of dynamic properties of supercon- 

10) is the time-dependent Ginzburg-Landau model^^ (TDGL). A 



ductors (see e.g. Refs. mi 

recent review- warns us of the fact that for gapped superconductors TDGL is not strictly 
valid even above T^, while below it becomes totally wrong; nevertheless TDGL remains 
popular because of its simplicity and its ability to reproduce observed phenomena. Kramer 
and Watts- Tobin^^ generalized TDGL so that it should be applicable to gapped supercon- 
ductors as long as there is local equilibrium, while still retaining some of the simplifying 
features of the TDGL formalism. 

The objectives of the present study are (i) to predict the influence of imperfections (i.e., 
nonuniformity and finite width) on fluctuation superconductivity in rings, (ii) to compare 
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our predictions with experimental results and thus check the range of applicability of TDGL, 
at least as a phenomenologic model, and (iii) to find how features of fluctuation supercon- 
ductivity (i.e., average current and lifetime of metastable states) are affected when TDGL 
is generalized to the Kramer-Watts- Tobin model. 

This article is organized as follows. In Sec. [ITl we describe our computational method, 
in Sec. IIIII we present general results that stem from our method and in Sec. IIVI we use 
our method to explain measured values. Most of our calculations aim at fitting the currents 
measured in Ref . [tI, but we also address experiments that measure the position of a vortex— 
and lifetimes of metastable states.- In each of these sections we distinguish between thin rings 
for which a one-dimensional treatment should suffice and wide rings for which corrections 
are necessary. In Sec. |V] we extend our analysis to the Kramer- Watts- Tobin model and 
the central features of our results are briefly summarized in Sec. I VII In Appendix |X] we 
describe how we dealt with some numerical difficulties and in Appendix [B] we speculate on 
the possibility of treating arrays of normal rings by means of the Ginzburg-Landau model. 



II. TDGL WITH THERMAL FLUCTUATIONS 



The Ginzburg-Landau model describes the behavior of a superconductor by means of an 
order parameter field ip and an electromagnetic potential A. For a thin wire and close to the 
onset of superconductivity, the magnetic field in the wire can be approximated by the applied 
field. In this case and in cgs-Kelvin gaussian units the model postulates the energy density 
a|?/'|2+/3|7/'|72+(nV2m)|(iV-27rA/<l>o)^/'P+A-j/c, with a = 12hkB{T-T^) /nmivF, ks the 
Boltzmann constant, T the temperature, Tc the critical temperature, m the mass of a Cooper 
pair, i the mean free path, vp the Fermi velocity, /3 = 87r(K^e/mc)^, n = 1.15 x — e 

the electron charge and j the total current density. The integral of the energy density over 
the volume of the sample will be denoted by G. All the experiments^'^'^ with which we 
will compare our results were performed on mesoscopic aluminum samples; we will therefore 
adopt vp = 1.5 X 10^ cm/sec, as appropriate for aluminum with a BCS coherence length 
1.6/im; this value of vp is smaller than what would be obtained within the free-electron 
model. We will assume that the samples are in the dirty limit. 

The characteristic temperature scale in our analysis will be Tlp = 7T^hvFi/2AkBL'^, where 
L is the perimeter of the ring. Tlp is the depression of the transition temperature that would 
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be obtained in the absence of fluctuations for a uniform profile and half-integer magnetic flux 
(in units of $o). Tlp equals the Thouless energy multiplied by tt/S/cb- For the conductivity 
of normal electrons we have taken the aluminum value a = 1.22 x 10^^£fi^^cm~^, but we 
have also checked other values. 

All the rings that we will consider will be thin and narrow in the scale of the coherence 
length i = h/ \j2m\a\ and the magnetic penetration depth k^. However, not all the rings 
will have a width which can be neglected when compared to its average radius; these rings 
will require a separate treatment. 



A. One-Dimensional Loops 

A numerical method for the implementation of TDGL with thermal fluctuations was 



presented in Ref. Il6| for the case of thin wires and loops that do not enclose a magnetic 
field. In this section we review this reference and point out the required modification when 
magnetic flux is enclosed. 

For computational purposes the loop is divided into segments of length L/N and the 
fields and the energy are discretized. A can be taken tangential to the wire and we denote 
by ipk and A}^ the values of ip and A in the A;**^ segment. We also define A^ = 2'KLAk/ N^q. 
From TDGL without fluctuations it follows that the derivatives of ipk and Ak with respect 
to time have the form 

dRe[iJk]/dt = -T^,kdG/dRe[iJk] , (1) 
with an analogous equation for Im[-?/'fc], and 

dAk/dt = -TA,kdG/dAk , (2) 

where F^^^ = Nmvpi/Sh'^Lwk and TA,k = ^e^L/Nh?aWk, with Wk the cross section of 
segment k. 

Equations ([1]) and ([2]) have the "canonic" form dx/dt = —TdG/dx. In this case the 
influence of fluctuations is obtained by adding Langevin terms to the right hand side of 
these equations, such that the integral of a Langevin term over a period of time r has 
gaussian distribution, zero average, and variance 2TkBTT. 

The electromagnetic potential can be gauged out of the expression for the energy, and 
therefore from Eq. ([1]), by means of the transformation ip{s) = U{s)il){s) with U{s) = 



exp[(27ri/$o) Jq A{s')ds']. The price of this transformation is that the gauge-invariant order 
parameter ip does not obey a "canonic" evolution equation. 

Reference ll6| adopts several normalizations that are appropriate for the fluctuation region 
T ^ Tc- The typical size of the order parameter is tp = {2mk%T^ / (3h?w'^Y^^ , where w is the 
average cross section, the typical coherence length is = {w^l/32TT^K'^kBTcy/^ and the time 
unit is f = Under free electron gas assumptions, the parameter 7 of Ref. Q can be 

expressed as 7 = 0.5(L/,^^)^. We then define the normalized quantities L = L(2mfcsTc)^/^/^, 
R = Ge'^L/amwvFi, 4'i3k = 4'k/'4', i^pk = 'ipk/i' and a' = a/f3ip'^. For the aluminum values 
and the units that we have adopted, = 0.679{w/K,^T^y^^, L = 6.72 x lO^LTc^^, R = 
4.62 X IQ-^^L/fw, a' = 8.55 x W^^jiT - T^)/iT^ and f = 2 x 10-^^1/^- With these 
normalizations, Eqs. ([T]) and become 

A^^ci'pk = -{T/f)[{a' + \tPpk\^)^Ppk + {N^p/L)\2wky^ 

[{Wk + Wk+l){ijpk - i^pk+l) + {Wk + Wk-l){'lplik - , (3) 

A^ac^fe = -{CAT/wk)[in/ekBT + {ip/L)\m[5k + , (4) 

where Ca = {wR/2Nf)[^i3L/ LY , I is the current around the ring, 5^ = N{wk-i + 
'Wk)4'*pi:-i'^l3k/u!, * denotes complex conjugation, 6n+i = ^i, wq = wn and ipi^o = 
exp(— 27rz$/<l>o)V'/3Af) where $ is the flux enclosed by the ring. Amac denotes the macro- 
scopic (i.e., it ignores fluctuations) increment during a short period of time r. 

The fluctuations in the increment of Ak have a variance ACA^/wk] the fluctuations of ipiBk 
are taken into account by adding to the real and to the imaginary part terms with variance 
Nrw^i^/fWkL and then modifying its phase in order to take into account the fluctuation in 
U. 

In this study we assume that nonuniformity enters the problem through the cross section 
of the ring; we expect that chemical nonuniformity would have a similar effect. 

Since the induced flux is negligible, if the applied flux enclosed by a ring remains constant, 
then Ak is also constant. This condition determines the current J; from Eq. (jlj) it leads 
to 

^ ^ ekeT ^ r]k/ r - {Ca^/b/ L)lm[Y^ w^^jh + 



h CaT.^. 



-1 

k 



where the sums are over the segments of the ring and rjk is the contribution of the Langevin 
term to Ak during the period of time r. 
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Quantities of interest are evaluated by explicit Euler-Mayurama steps following Eqs. ([3]) 
and (jl]). We start from an arbitrary set of values for the order parameter, the initial steps are 
intended for relaxation to some typical state, and subsequent steps are used for statistical 
averaging. 

The only modification that is required in order to apply the procedure of Ref. 
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IS m 



the boundary condition for ip, which in the case that flux $ is enclosed becomes ipi^) 
exp(27ri$/$o)^(0). 



B. Wide Rings 

Although theory is "cleaner" for ID loops, in real life this idealization might not be 
justified. Moreover, wide rings carry larger currents and may provide us with reliable mea- 
surements in fluctuation regions of the $ — T plane, where thin rings have a low signal to 
noise ratio. It will therefore prove useful to develop a simple model for rings that are not 
ideally ID. 

We consider superconducting rings with a width that is very small compared with the 
lengths over which the superconducting variables change significantly (coherence length, 
effective magnetic penetration depth), but is not sufficiently small compared with the radius. 
As a consequence, the fluxes enclosed by the inner and by the outer boundaries will not be the 
same. Nevertheless, we will build a quasi-lD model in which the superconducting variables 
will be represented by quantities that do not depend on the radial coordinate. 

Let us denote by R the average radius and the width by D{6); the thickness is w{6)/D{6). 
The ring occupies the region R — D/2 < r < R + D/2 and is immersed in a uniform 
perpendicular magnetic field B. At this stage we neglect the induced magnetic field and B 
will be taken as constant. 

With the purpose of averaging over r, we separate the electromagnetic vector potential 
into a contribution of the applied magnetic field and a contribution a due to thermal fluc- 
tuations, i.e. we write as A = {Br/2 + a)6. In order to have a ID model we assume that ijj 
and the voltage are functions of 6 and time only, which implies 

R R ^ 

a = an—, V^" = i^R—O , (6) 

where an is the value of a at r = and ip'j^ is the derivative of ip with respect to the arclength 
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s = R6, taken along the circle r = R. We will also approximate the current density by the 
form j = jji{R/r)6, with jn independent of r and proportional to l/w{6). 

We now integrate the energy density over R — D/2 < r < R + D/2 and neglect terms 
of order higher than {D/Ry. In this approximation the volume averages are ((r/i?)") = 
1 + n{n + l)(D/i?)^/24 and the free energy becomes 

G = ^wds{a\tlj\^ + /3\^\y2+{ff/2m)[{27r/^of\ijWB\AR^ + D^)/16 
+aRBR + alfD) + WfD 

+ {2m/%)iBR/2 + aMii^^P'R - + iBR/2 + fDaR)jR/c} , (7) 

where /d = 1 + D^UR^ + 0((D/i?)^). 

We can now introduce a new variable, A = (27r/$o)(-Bi?/2//) + and rewrite G in the 
form 

G = ^ wds{[a + {eBDf/{Qmc^)]\iPf + 131^1^2 + 

ifnhy2m)[\i;'^\' + i^l' + ^M^^'*r ' + fD%AjR/27rc} . (8) 



Discrete variables can now be defined as in Ref. 
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: ijjk as the average of in cell k and Ak 



as the integral of A along the length L/N of this cell. It can be checked that these variables 



behave as "canonic" in the sense of Ref. 
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Integrating over r the TDGL expression for 



dip/dt, under natural averaging assumptions we recover Eq. ([T]), with G given by Eq. ([8]). 
However, using {^Q/27i)dA/dt = daR/dt = (dA/dt) = —c{jN)/o', where Jn is the normal 
current density, we obtain that the coefficient F^i ^ is smaller by a factor fu than the expres- 
sion obtained for a ID ring. The same simplifications that were obtained in Ref. for ID 
loops can be achieved here by defining the gauge-invariant order parameter ip^s) = U{s)ip{s) 
with U{s) = exp (^i A{s')ds'^ . 

By detailed comparison we conclude that a ring with non-negligible width can be treated 
by means of Eqs. ([3]) and (jl]) and the fiuctuations described under them, provided that we 
include the following corrections: (i) a has to be replaced with a + {eB Dy / [Gmc^) , which 
means adding (l/3)($D^^/$o-R^)^ to a'; (ii) the term proportional to ^'^ in Eq. ([3]) has to 
be multiplied by fo and (iii) the variance of Ak has to be divided by a factor Jd. 

As long as we neglect the self inductance, the integral of the electric field (and hence of 
daR/dt) around the ring has to vanish; in this case the instantaneous current is given by 
Eq. ([5]). If there is non- negligible infiuence of the ring self inductance, which we denote by 
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Ls, the total flux through the ring will not be the applied flux but rather $ = ^x + cLgl. 
This has two consequences. The flrst is that f A ■ ds varies with I and Eq. ([5]) generalizes 
to 

i2eLs/hT) + (hCA/ekBT) Z w^' ' ^ ^ 

where /prev(^) = I{t — t). For the range of parameters that we considered, this reflnement 
turned out to have no noticeable effect. The second consequence is that when comparing 
our results with experiments, we should present the results as functions of = $ — cLsI, 
since usually this is the controlled quantity. For the parameters we considered, |<l>a; — $1 is 
small but noticeable. 



III. REPRESENTATIVE RESULTS 
A. One-Dimensional Loops 

Our formalism enables us to study samples with cross section given by any periodic 
function w{6) of the angular polar coordinate by 6. Therefore, our flrst question is what 
nonuniformity proflles would it be interesting to study. Since any periodic function can be 
written as a Fourier series w{9) = w{l + ^ (3j cos j9 + ^ 7^ sin j6'), another way to pose the 
question is what Fourier coefficients (3j and 7^ would it be interesting to consider. 

We will restrict ourselves to rings with mirror symmetry w{—9) = w{6), so that 7j = 0, 
and will have to decide on the choice of the /3j's. We may be guided by our results for rings at 
T ^ Tc and small deviations from uniformity, in the absence of fluctuations.— »^>^ In those 
studies, which were corroborated by other groups, ^^'^^ it was found that nonuniformity can 
qualitatively modify the Little-Parks phase diagram and there is a range of temperatures for 
which there is a continuous passage between consecutive fluxoid states. At the lower end of 
this range, there is a critical point P2 such that below P2 the passage between fluxoid states is 
hysteretic and such that at P2 the derivative of the current with respect to the flux diverges. 
The important aspect of our previous results for the question at hand is that the salient 
features of the phase diagram, such as the position of P2, depend only on the eccentricity 
Pi, whereas higher harmonics have no influence at the leading order. It is therefore natural 
to consider a shape in which only the flrst harmonic is present, i.e. w{6) = w{l + Pi cos6). 
This proflle will be called "sinusoidal." 
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FIG. 1: Cross section profiles considered in this study. 



In order to gain some understanding of how other profiles behave, we consider cross 
sections with discontinuous functions w{9), since in this case the coefficients have the 
slowest decrease with j, providing a complementary situation to the sinusoidal profile. The 
simplest discontinuous profile is that of a ring with a constriction. This situation might also 
be of experimental interest, since the constriction might represent a Josephson junction or 
a defect. Finally, in order to have an opposite case to the sinusoidal profile, we consider 
the case of rings with an additional mirror symmetry w{'k — 9) = w{6), since in this case 
Pi = 0. A profile like this will be called "symmetric." Figure [1] shows the nonuniformity 
profiles which we have considered. 

Figure [2] shows the current J as a function of the magnetic flux $ for fixed temperatures. 
The choice of geometric and material parameters was inspired by those of the samples for 
which data are available.- The upper panel is for T = Tc—l.GiTi^p and all the curves represent 
sinusoidal cross section profiles. We see that nonuniformity has two qualitative effects: one 
of them is overall decrease of the current and the second is a shift of the maximum current 
to lower fluxes, so that the decay of J as $ — ?■ $o/2 becomes more gradual. In contrast 
to the well known behavior in the absence of thermal fluctuations, even for uniform cross 
sections and T < Tc — Tlp, the average current vanishes at $ = $o/2 and the passage 
between consecutive winding numbers is smooth. 

The lower panel in Fig. [2] is for T = Tc — 2.67Tlp. For a uniform sample at this tempera- 
ture, the average current does not go to zero for $ — t- $o/2. According to a blind application 
of statistical mechanics, this average should vanish, since the energies for clockwise or anti- 
clockwise currents are equal; however, there is an energy barrier between the two fiuxoid 
states, so that in real experiments or simulations only one of them is probed. For all the 
nonuniform samples shown in the graph, the current does go to zero as $ — t- $o/2, meaning 
that either there is no energy barrier, as is the case for temperature above that of P2, or the 
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barrier is not large in comparison to thermal fluctuations. 

The profiles with a constriction in Fig. |2]had a segment of length 0.8L with large cross 
section and a constriction of length 0.2L. We see that the /($) curves for these profiles 
practically coincide with those of appropriate sinusoidal profiles. (Only one line was drawn 
in these cases for both sets of symbols.) From here we may adopt the working assumption 
that for any reasonable profile there will be an equivalent sinusoidal profile, so that by studing 
sinusoidal profiles we may expect to obtain most of the interesting information. Defining 
the eccentricity as the first harmonic of the cross section divided by its average, we obtain 
that the constricted profile equivalent to Pi = 0.35 has eccentricity 0.22 and the sample 
equivalent to f3i = 0.8 has eccentricity 0.4. It follows that the /($) curve is influenced most 
significanly by the eccentricity, but not only by it. As an extreme example, we considered 
a symmetric sample for which the eccentricity is zero. This sample consists of two wide 
segments, each of length 0.2L, connected by constrictions of length 0.3L. In spite of the 
very large ratio between the cross sections, which manifests itself in a strong inhibition of 
the average current, the position of the maximum remains quite close to $ = $o/2. As a 
general observation we might say that the inluence of nonuniformity on /($) is stronger for 
lower temperatures and close to $o/2. 

Figure [3] shows the ratio of the current at a given fixed flux to the current at flux $ = 
O.l'I'O) as a function of temperature. Again, we see that the influence of nonuniformity 
increases as $ approaches $o/2 and, for the parameters used in this graph, is noticeable for 
T < Te - I.ITlp. 

B. Wide Rings 

We consider first the case in which the cross section is uniform along the ring. We denote 
by Id the current around the ring for non-negligible width, and by Jq the current in the 
case D— 7-0,^//^— T-oo (the cross section is the same in both cases). In all the cases we 
considered, we found < /d/Jq < 1, i.e., the finite width inhibits the current, but does not 
wipe it out. The inset in Fig. H] compares Id and Iq for a typical situation. 

In order to make more quantitative statements, we classify our results into those outside 
and within the fluctuation regions. In the absence of fluctuations, our quasi- ID model leads 
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FIG. 2: Average current as a function of the magnetic flux for fixed temperatures and various 
cross section profiles. In all cases we took L = 27r x lO^^cm, iv = 6.38 x lO^^^cm^, Tc = 1.251K, 
i = 2.66 X lO^^cm, 10 computational cells, 2 x 10® averaging steps, 6 x 10^ relaxation steps, each 
step of 10^® time units. Meaning of the symbols: ♦ uniform cross section; sinusoidal with 
Pi = 0.35; ■ sinusoidal with Pi = 0.8; Q profile with constriction; A symmetric profile. For the 
profiles with constriction the cross sections of the narrow and the wide parts are in the ratios 1:2.14 
and 1:7.97; for the symmetric profile they are in the ratio 1:27.3. The curves are guides for the eye. 

to vanishing order parameter when 

(l/3)(D<l>/i?<l>o)' + /d($/$o - = (T, - r)/4TLP , (10) 

where n is the winding number. Let us denote by $„_ the smaller (larger) root of 

Eq. ( fTOi) : the regions < $ < where the sample would be normal in the absence 

of fluctuations, will be called "fluctuation regions." 

Figure m shows the difference Id — Iq outside the fluctuation regions for two temperatures 
and two widths, normalized by the factor {D/Ry. The four sets of points lie in a nearly 
universal curve, indicating that for the tested range Id — h is proportional to and its 
value is independent of Tc — T. 

Figure [5] shows the current within fluctuation regions for the range |$ — $n±| ^ 0.06$o- 
If the flux is measured from $^±5 Id{^) is almost independent of D (but not of the tem- 
perature) for D/R < 0.7 (respectively 0.3, 0.2) in the case ^n± = '^'o+ (respectively 
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FIG. 3: Current ratios /(0.4$o)//(0.1$o), /(0.45$o)//(0.1«>o) and /(0.4875$o)//(0.1$o) as func- 
tions of the temperature, for nonuniformities given by /3i =0, 0.35 and 0.8. The parameters and 
symbol meanings are the same as in Fig. [2l The circles are for a sample with a constriction of 
length 0.2L and cross section equal to 1/2.14 of that of the rest of the sample; they are shown 
for I{0.4:5^q)/I(0.1^q) only. The calculated symbols have been joined by straight lines. For visi- 
bility, the curves for /(0.45<l>o)/-^(0.1*&o) (i"ed online) have been raised by 1.5 units and those for 
/(0.4<I>o)//(0.1<I'o) (black) have been raised by 3 units. 



Let us now consider nonuniform cross sections. We have to deal separately with the case 
in which nonuniformity is due to the thickness, so that D is constant, and the case in which 
nonuniformity is due to the width, so that D is a function of 6. For the parameters we 
considered, there was no significant difference; the results we report in the following are for 
uniform D. Figure [H] compares the shapes of the current-flux curve for a uniform and a 
nonuniform sample. As in the case of ID loops, nonuniformity brings about a more gradual 
decay at $ ~ 0.5<l>o- 

We would finally like to know whether the influence of eccentricity is enhanced or 
weakened by the width of the ring. Denoting by l£){Pi,(f)) the current for a ring of 
width D with eccentricity f3i at flux 4>^o, we may regard the ratio 0.45)//d(/3i, 0.1) 

as a representative index for the shape of the function //)(/3i,0), and the difference 
/£)(0.7, 0.45)//£)(0.7, 0.1) — /^(O, 0.45)//i3(0, 0.1) as a representative index of the influence 
of eccentricity on the shape of Id{Pi,4>)- We see from the inset in Fig. [6] that for wider 
samples the influence of nonuniformity is smaller. 

In our calculations we used several values for the conductivity, with a ratio ~ 3 x 10^ 
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FIG. 4: Reduction of the current as a function of the magnetic flux, outside the fluctuation regions. 
Calculations are for a sample with R = 3.5 x lO^^cm, w = 7.83 x lO^^'^cm^, Tc = 1.24K, i = 
3.02 X 10~^cm; computational parameters as in Fig. [2j Symbols: □ T = Tc — 0.67Tlp, D/R = 1/3; 
X (red online) T = - 0.67rLP, D/R = 2/3; (blue online) T = - 1.17rLP, D/R = 1/3; 4- 
T = Tc- 1.17rLP, D/R = 2/3. The value for T = - 0.67Ti^p, D/R = 2/3 and ^> = 0.65$o lies 
within the fluctuation region 0.38 < ^/^o < 0.70 and has been included for comparison. Inset: 
currents for ID and quasi-lD samples as functions of the flux; o T = Tc — 0.67Tlp, D = 0; x (red 
online) T = - 0.67rLP, D/R = 2/3. 
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FIG. 5: Current within fluctuation regions in the vicinity of the region border, ^n±, for two fixed 
temperatures. For visibility, the values for T = — 0.33Tlp have been shifted 0.02 units to 
the right. Symbols: □ ^>„± = ^>o+, D/R = 1/3; x ^>„,± = $o+, D/R = 2/3; + $„± = 
L)/i? = 1/3; * <I)„_I_ = <I)^^, D/i? = 1/5; o L) = 0, arbitrary <I>„±. Other parameters as in Fig. [H 
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FIG. 6: Current as a function of the flux for samples with width D = R/3 at T = Tc — 1.17Tlp. 
The curves are guides to the eye. The rhombs (joined by a blue online curve) are for a uniform 
sample, and the squares are for w{0) = ■w{l + 0.7 cos 6). The other parameters are as in Fig. SI In 
order to compare the shapes, the results for f3i = 0.7 have been multiplied by 1.37. The inset shows 
an index that represents the influence of eccentricity on the shape of the current, as a function of 
the width of the sample, "ioe" stands for /z)(0.7, 0.45)/Id(0.7, 0.1) - /d(0, 0.45)//d(0, 0.1). 

between the largest and the smallest. The value of a had no appreciable influence on the 
average current. This can be understood, since the average current is an equilibrium quantity 
and a only affects the dynamics. 

IV. COMPARISON WITH EXPERIMENTS 
A. One-Dimensional Loops 

We have used our method to analyze the 15 samples studied in Ref. Q.— Samples for 
which D < R/A were treated as "one-dimensional." We modeled each sample by a ring with 
sinusoidal cross section profile. For a given temperature, our first task was the generation of 
a representative curve /($) for the experimental data. As the first step we subtracted the 
linear background; in some cases, we assumed that for a ID sample the maximum value of 
|/| in the range < |$| < 0.5$o has to be the same as for 0.5$o < 1"^! < "^o- We then shifted 
the origin of the applied field to inforce f^^^^^^^I^^) cos{27r^/^o)d^ = 0, and finally fitted 
the data in the range — 0.5$o ^ ^ ^ 0.5<l>o to a high order odd polynomial that vanishes at 
$ = ±0.5<l>o. 
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The second task was the choice of parameters Tc, £, /3i and the mutual inductance between 
the sample and the scanning SQUID. Tc was taken from Ref. [T). i was chosen by assuming 
a uniform ring, considering the highest temperature for which measurements are reasonably 
reproducible and periodic, and then requiring that our calulated values reproduce the exper- 
imental smoothed ratio /(<l>i)//(<l'2) where, typically, $i ~ 0.1$o and $2 ~ O.SS'l'o- In order 
to fix /?!, we again require agreement with the experimental ratio /($i)//(<l'2), but this time 
we consider the lowest temperature for which there is no hysteresis and $2 ~ 0.47$o- After 
/3i was fixed, we re-evaluated i, using the obtained value of f3i rather than /3i = 0; we found 
that this second iteration did not modify i significantly. Finally, the mutual inductance was 
fixed by fitting the calulated points to the experimental /($) for the lowest temperature 
and in the range < $ < 0.4$o- The mutual inductances we found are typically (10 ± 5)% 
smaller than those reported in Ref. Q. This small discrepancy can be attributed to the 
experimental uncertainty in the sample-SQUID distance. 

For most of the samples, the best fit was obtained when a uniform cross section was 
assumed. However, there were a few samples for which a better fit was obtained for /3i 7^ 0, 
and two of them are presented in Fig. [71 each for a temperature above or close to Tc— Tlp and 
another temperature well below — Tlp. The temperature of P2 is ~ T^ — (1 + 2/3i)Tlp,— 
so that we may expect nonuniformity to extend the range over which there is no hysteresis 
at $ = $0/2 by an amount of ~ 2/3iTlp. In Fig. [3, the lines marked with a temperature 
are the experimental representative curves for /($), the points are calculated values using 
our method, and the lowest line is an interpolation for values calculated assuming a uniform 
cross section. Sample 13 is the one with largest deviation from uniformity and, indeed, SEM 
inspection reveals that this ring has imperfections. In general our calculated points are in 
good agreement with the experimental lines. For the lower temperatures shown in Fig. [71 
the theoretical curves are more rounded than the experimental curves; this is probably due 
to our model assumption that the cross section profile is sinusoidal. 
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FIG. 7: Current as a function of the magnetic flux for two samples that appear to be nonuniform. 
The lines marked by temperatures are experimental and the points were calculated. The lowest line 
(blue online) was calculated for /3i = 0. Parameters for sample 7 taken from Ref. lZ|: L = 3.14/im, 
w = 3.8 X lO^^^cm^, Tc = 1.264K; parameters for sample 7 fitted here: £ = 20.5nm, /3i = 0.13 and 
mutual inductance Ms-r = 0.079/i$o/nA. For the line /3i = we took Ms~r = 0.078/x$o/nA. 
For sample 13, L = 6.28//m, tt) = 6.4 x lO^^^cm^, = 1.251K, I = 26.6nm, j3i = 0.35, Ms-r = 
0.30^<l>o/nA; for /3i = 0, Ms-r = 0.27/x<I>o/iiA. The kink in the blue line for sample 13 is an 
artifact of the algorithm that interpolates between our calculated points. 

B. Wide Rings 

1. Current-flux curves 

In this section we deal with the samples of Ref. Q that have D > R/4. As in the case 
of ID rings, we generated representative curves /($) for the experimental data, but slightly 
different strategies were adopted: the overall slope of the data was not subtracted, the shift 
in the origin to render /($) an odd function was fixed "by eye," and different polynomia 
were used to fit the data in different ranges. 

As discussed in Sec. IIIBl the contribution of the induced magnetic field is not quite 
negligible. Since the field generated by self induction is not uniform, and since additional 
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magnetic induction mechanisms may be present besides self induction, we initially attempted 
to regard Lg as an adjustable parameter. The results we obtained were of the order of the 
tabulated values,— but systematically smaller; in some cases we obtained < 0. In view 
of this behavior, and in order to reduce the number of parameters, we decided to neglect 
self induction altogether. This omission may influence the effective value of 

Figure [8] shows our results for sample 5. This sample had a width/radius ratio D/R = 
0.39. Note that for this sample the vOR method could not be applied.- We chose the value 
of i from the data for T = Tc — 0.17Tlp, the value of f3i was adjusted to reproduce the 
experimental ratio /(0.507$o)/-^(0.1$o) at T = Tc — 1.56Tlp and the mutual inductance was 
adjusted aX T = — 1.17Tlp. We also used the data at T = Tc — 1.56Tlp to refine the 
calibration of the applied flux. Our calculated results in the hysteresis region (inset) are 
in surprising agreement with the experiment. However, the calculated decays of metastable 
states are much sharper than the observed decays; the lack of sharpness in the observed 
decay is probably due to the fact that the experiment was not perfectly static: the range of 
flux involved in the passage between fiuxoid states was swept during a lapse of time of the 
order of 10~^ sec, which is comparable with the filtering time. 

Figure [9] shows our results for sample 12, which had a width/radius ratio D/R = 0.35. 
Among the samples that we consider as wide, this is the one for which the largest flux was 
reached in the experiment. The value of i was adjusted using the data for T = T^ — 0.42Tlp 
and (3i and the mutual inductance were adjusted using the data for T = Tc — 1.42Tlp. We 
see that there is fair agreement between our model and the experimental results for $ < $0; 
but for larger fluxes there is a considerable deviation. The experimental currents vanish for 
$ ~ 0.98$05 whereas the calculated currents vanish almost exactly at $ = $o- 

The fluctuation region around $ = 0.5$o for sample 12 is shown enlarged in Fig[TOl In 
the absence of fluctuations, the current would vanish for 0.32 < $/$o ^ 0.71. The lines 
that show the current that would be present without fluctuations were obtained with the 
same code as all the other results, but the Langevin terms were divided by a factor of 100. 
Note that all the adjustable parameters of our model were fixed outside this region. 

Our fitted values of £ vary from sample to sample. The largest value equals almost twice 
the smallest value, but they are all (including thin and wide samples) within 15% of the 
empirical expression 1/i = 1.5 x lO^cm"^ + 2/D + 14.5/L. Our values of i are larger than 
those of Ref. Q; the difference is mainly due to the adopted value of vp- 



18 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 



$/<Do 

FIG. 8: Current as a function of the magnetic flux for a wide sample. The lines are smoothed 
experimental data and the dots were calculated. The larger currents are for T = Tc — 1.17Tlp and 
the smaller currents (blue online) for T = Tc — 0.67Tlp. Parameters taken from Ref. 7: L = 2.2/im, 
w = 7.83 X lO^i^cm^, D = 1.35 x lO^^cm, = 1.240K; parameters fitted here: I = 30.2nm, 
/3i = 0.07 and mutual inductance Ms-r = 0.035//$o/nA. Inset: hysteresis region; the curves 
(actually, dense set of points, red online) are experimental data for T = Tc — 1.56Tlp and the black 
lines are calculated. 
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FIG. 9: Current as a function of the magnetic flux for another wide sample. The larger currents 
are for T = Tc — 1.42Tlp and the smaller currents (blue online) for T = Tc — 0.42Tlp. For visibility, 
the currents in the case T = Tc — 0.42Tlp have been multiplied by the factor 3. Parameters taken 
from Ref. ^: L = 3.14/xm, id = 1.01 x lO^i^cm^, D = 1.75 x IQ-^cm, Tc = 1.244K; parameters 
fitted here: i = 30.9nm, /3i = and mutual inductance Ms-r = 0.082/x<I>o/nA. 
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FIG. 10: Blowup of the fluctuation region around ^> = 0.5^>o for T = Tc - 0.42rLP in Fig. El The 
obhque Unes of straight appearance at the left and at the right describe the currents that would be 
obtained without thermal fluctuations. The curve describes the smoothed experimental results and 
the symbols were calculated. • TDGL; x KWT (Sec. lY]), Tph = lO'^ sec; + KWT, Tph = IQ-^ sec. 

2. Flux-induced vortices 

A recent experiment in a doubly connected asymmetric disk— found that, for the 
examined sample, the order parameter practically vanishes in the narrow part of the sample 
when passing between the fluxoid state and the fluxoid state 1; when passing between 2 
and 3 the order parameter practically vanishes in the wide part, and no local vanishing is 
observed in the passage between 1 and 2. 

We consider now a sample with geometric and material parameters similar to those 
reported in the experiment and check whether our theoretical method reproduces the exper- 
imental behavior. Our results are shown in Table HI We denote by ipe the order parameter at 
angle 6, where ^ = is the angle at which the sample is widest; (■ ■ ■ ) is average over time; for 
each temperature and flux we use the normalization constant C = 0.5((|'?/^^/2p) + (|'?/'_,r/2p)); 
stdevg stands for (dV'el^) ~ (I'^eP)^)"'^^^- Indeed we see that for $ = 0.5$o and for a tem- 
perature within the appropriate experimental range (IV'P) is particularly small for 6 = ir, 
whereas for $ = 2.5$o the small value is obtained for 6 = 0. In the table we show the first 
two moments only; higher moments indicate that in these two situations {ipl'^ has exponential 
distribution. 

Moreover, for $ = 1.5$0) (l^eP) is not small for any 6, in agreement with the experimental 
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TABLE I: Angular dependence of the characteristic size of the order parameter for a sample with 
radius 290 nm, D{e) = 168(1 + 0.625 cos 6^) nm, w = 5.04 x lO^^^cm^, = 1.36K and £ = 15.6 nm. 



C = 0.5((|V'7r/2p) + (l^-7r/2p))- Appropriate temperatures were taken fron Ref. 
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FIG. 11: Average value of IV'P as a function of the position in the sample for <I> = 1.5<I>o- The 
sample is widest at 6 = and narrowest at 9 = ir. The sample properties are described in Table [H 
□ x) T = 1.19 K (T = 1.2 K, T = 1.21 K), 12 computational cells. Typically, we took 3 x 10^ 
averaging steps, 6 x lO'' relaxation steps, each step of 3 x 10^^ time units. The lines are guides for 
the eye. Note that the normalization C depends on the temperature. 

observation. Figure [TT] shows that for this flux and temperature (iV'eP) does not depend 
very strongly on 6. However, we note that stdevg (respectively stdev^r) is almost as large as 
(I'^/'oP) ((IV'ttP))) suggesting that there are frequent transitions in both directions between 
the fluxoid states 1 and 2. Since a continuous change in winding number must involve a 
place where the order parameter vanishes, we may anticipate that small averages of [t/'P will 
be obtained if we average over transition steps only. 

For this purpose we need a criterion to decide what is meant by a "transition step." We 
start by noting that in the fluxoid state 2 a positive current flows around the sample for 
$ ^ 1.5<l'0) whereas in the state 1 the current is negative, so that in a transition between the 
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FIG. 12: Probability per unit angle for the location of the vortex when the fluxoid state of the 
sample switches between 1 and 2. * T = 1.22 K; other parameters and symbols as in Fig. [TTJ 

two states the current should change sign. However, there may be many small steps in which 
the sign changes back and forth; we will not count each of these changes as a transition, but 
only the first one after the sample has been in a "typical" state. Let us denote by /+ (/_) 
the average of the current over those steps in which it is positive (negative). We found that 
|/_| ^ /_|_. We regarded a step as a "transition up" (down) if it is the first one for which 
the current is positive (negative) after having been smaller than —1^/8 (larger than /+/8). 
When a transition step is detected, we identify the cell in which lipl"^ has its minimum value 
as the place where a phase slip occurs. 

Figure [12] shows the probability density for the location of the phase slip in a fluxoid tran- 
sition as a function of the angle 6, for $ = 1.5$o and for several temperatures close to the 
onset of superconductivity. There is no appreciable difference between the probability dis- 
tributions for transitions up or transitions down, and therefore we present the distributions 
for transitions in both directions. The probability distributions do not depend appreciably 
on the criterion for "typical state"; for instance, using the threshold /+/4 rather than /+/8 
leads to practically the same distribution. The probability distributions slightly depend on 
the value of the conductivity; the effect of raising the value of a is similar to that of lowering 
the temperature. The results in Fig. [T2] were obtained for a = 1.22 x lO^^iQ^^cm^'^. 

As expected, transition steps occur when there is at least one cell where is small. If 
we denote by A the event that there is a cell where < C/20 and by B the event that 
there is a transition step, then, for T = 1.2 K and $ = l.S'l'O) the probability of having 
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simultaneously A and B is 1.9 times larger than the product of the probabilities of A and 
B. Transition steps are clustered: denoting by ttr the time interval between consecutive 
transition steps and by overline the average over intervals, the relative standard deviation 
is (t^j. — ttr^)^^^/^tr = 3.7. This result is consistent with a scenario in which most of the time 
the order parameter is too large to permit transitions and, every time it becomes suficiently 
small, several transitions occur. 



3. Life-time of metastable states 

Zhang and Priced performed direct measurements of the lifetime of metastable states, as 
a function of the temperature and the flux. They prepared states with winding number 
that enclosed flux larger than 0.5<l>o and waited until a change in the magnetic susceptibility 
was detected. The elapsed time was measured and the same procedure was repeated several 
times, until a significant average was obtained. 

There are many reasons not to expect that our method will be able to reproduce these 
experimental values of the lifetime. First, TDGL is not expected to be valid at the considered 
temperatures; second, the width/radius ratio for their sample was 0.54, so that the quasi- 
ID description seems inappropriate; third, the direct measurements could be performed 
for lifetimes larger than 1 sec, whereas the lifetimes that we could practically study were 
smaller than 10"'^ sec; fourth, since the fluctuations above in Ref. g] were much larger than 
predicted by theory, we may suspect that some non-thermal perturbation was present in the 
experiment and that this perturbation (rather than thermal activation) could induce the 
decays. In spite of all these difficulties, having an estimate of the influence of nonuniformity 
on the lifetime of a metastable state is at least of academic interest. 

Let us denote by 0_4 the value of for which the lifetime of the state with winding 

number is 10~^ sec. The values of 0_4, for temperatures at which measurements were 
performed and are not too far from Tc, are shown in Table [Tll The experimental value of 
0_4 was determined by extrapolating the exponential dependence found in Ref. ^ For com- 
parison, 0_4 was also calculated with our method, using the physical parameters reported 
as the best fit in Ref. and two cross section profiles; in one case the cross section was 
taken as uniform and in the other case we assumed uniform width, whereas the thickness 
had sinusoidal dependence, with /3i = 0.25, which is not far from the value 0.15 used for the 
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TABLE II: Flux for which the hfetime of the metastable state is 10~^ sec. The values of (/)_4 in 
the 2^^^^ column are experimental and in the 3^'^ and 4*^^ column were calculated using TDGL; 
the last two columns are for the Kramer- Watts- Tobin model (Sec. IV]). Physical parameters: 



radius=950nm, D = 


515 nm, w = 1.03 x lO^^^cm^, Tc = 


1.266K, a = 3.5 


;8 X lO^O-^cm-i, 


e = 31.8 nm, Tph = 10" 


sec. = 5; the time step was 3 x 10" 


time units. 




T(K) ,/._4 (exp.) , 


^TDGL ^ 0) ^TDGL ^ o_25) 


<ATM/3i = 0) < 


^KWT ^ 0.25) 


1.192 0.63 


0.72 0.62 


0.79 


0.62 


1.137 0.79 


0.88 0.73 







same sample in Ref. in order to mimic the shape of the phase diagram. The initial values 
of the order parameter were those that minimize the free energy for winding number 0; evo- 
lution followed governed by TDGL with fluctuations, and we decided that the initial state 
had decayed when a positive current was reached; after a decay, the process was repeated, 
during a total of 3 x 10^ time units. 

Although in Sec. IIVB II we ignored self inductance, in the present case flux sensitivity is 
very large and we therefore decided to take self inductance into account. Self inductance was 

n 

taken into account also in Ref. la, but the influence they found had opposite sign than what 
we find here: in our (resp. their) case self inductance increases (resp. reduces) metastability. 
In our case the main effect of self inductance is due to the difference between $a, and $; in 
their case the self inductance is used in a model for estimating the height of the potential 
barrier. 

Calculated lifetimes divided by their averages have typically Poissonic distributions. 
While bearing in mind the reservations raised at the begining of this section. Table IXTl 
suggests that nonuniformity could be a reason for the low lifetimes found in Ref. id. 

V. BEYOND TDGL 

A. The Kramer Watts-Tobin model 

Kramer and Watts- Tobin^^ (KWT) extended TDGL so as to render it applicable to 
gapped superconductors and valid as long as there is local equilibrium. The KWT model 
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has been successfully used to describe the current- volt age characteristic of thin wires.— — 
The energy functional is the same as in the Ginzburg-Landau model; in a gauge such that 
the electrochemical potential is uniform, Eq. ([T]) generalizes to 

1 

and similarly for Im[^/'fc]. Here K = 2A7 n'^e'^kBTcivFT^i^/hmc^, where Tph is the electron- 
phonon inelastic scattering time. In the limit K ^ 0, TDGL is recovered. Aluminum is 
challenging since it has a particularly large Tph, i.e., Xph = 10~^sec. As a consequence, the 
expected correction to TDGL should be large. Also, the range of validity of KWT should 
be small, Tc — T < h/kBT^h ~ 1 niK. 

Equation (fTTj) can be brought to canonic form by writing ip^ = \iJk\ exp(ixfc)- We obtain 

d\i,k\/dt = -h\.^\{i,k)T^^,kdG/d\i,k\ , (12) 
dXk/dt = -h^{i)k)^i,,kdG/dxk , (13) 
with h\^\{iljk) = {1 + K\iJk\^)-^^^ and h^{tPk) = 1/ {h\^\\iJk\^). 



d 

Jt ^ 



K_dM 
2 dt 



dG 



ipjk 



(11) 



B. Appropriate Langevin terms 



Following the reasoning of Ref. |16|, for a period of time r over which lip^l does not change 
significantly, the fluctuating parts of the changes of \tpk \ and of Xk would be expected to have 
variances {r]'^^^) = 2h\^\{il)k)T^^kkBTT and (r/^) = 2T^^kkBTT /h\.^\{il)k)\iJk? . However, we 
have pointed out elsewhere^ that in the case of Eq. (fT2l) this is not the full story. Due to the 
I 'j/'fc I -dependence of h\^\ and of the Jacobian W = d{Re[ipk],^^['ipk]) / ddiJkl, Xk) , fluctuations 
of lipkl do not have zero average but rather {ri\^\) = [dlog{h\^\W)/d\ipk\]h\^\T^^kkBTT. Using 
the expressions for and W we obtain 



\'ipk\(t + t) - \ipk\(t) = h\^\T^^k 



h^^^kBT dG 



where fj\^\ has gaussian distribution, (f^i^i) = and (^^^|) = '^h\^\{jpk)V^^kkBTT. 



(14) 



C. Results and comparison with experiments 

We have limited our study to two experiments. The first is the case of persistent current 
in a wide ring in the fluctuation region, shown in Fig. [101 The oblique crosses are for 
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FIG. 13: Increment of the flux for which the lifetime of a metastable state is 10"^ sec for T = 
1.192 K, /3i = (blue onUne) and /3i = 0.25 (violet online), for < Tph < 2.5 x 10"^° sec. The 
other parameters are as in Table HIl The dots were calculated and the lines are empiric fits. 

Tph = 10^^ sec and the upright crosses for Tph = 10^*^ sec; the other parameters are as 
for TDGL. We see that even for large values of Tph the currents obtained coincide within 
statistical uncertainty with those obtained with TDGL. 

The second experiment we studied was that of the lifetime of metastable states, considered 
in Sec. IIVB 3[ Figure [13] shows our results for the flux at which the lifetime is 10~^ sec. The 
results are reasonably fltted by expressions of the form 0_5(rph) = 0-5(0) + a[(6^ + r'^^^Y^^ — 
a/6], where a and h are adjustable constants. Using a flt of this form for 0_4 and extrapolating 
to Tph = 10~* sec leads to the results shown in Table HTl (The increment of 0_4 due to KWT 
for /3i = 0.25 is less than 0.01). 

VI. DISCUSSION 

We have evaluated the influence of thermal fluctuations on persistent currents (i.e., av- 
erage current), the position of the flux- induced vortices, and lifetime of metastable fluxoid 
states for aluminum rings that enclose magnetic flux, in a range down to ~ 10^ mK below the 
critical temperature, using TDGL with the addition of Langevin terms. In the cases that we 
considered to be most interesting, the evaluation was also performed using the KWT model. 
The considered rings were not "ideal," in the sense that they were either wide, so that the 
"enclosed" flux was not sharply deflned, or nonuniform, so that expansion of the order pa- 
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rameter into a Fourier series does not simplify the problem. In all cases we obtained at least 
qualitative agreement between TDGL and the experimental results, and in many cases we 
obtained quantitative agreement, in spite of the fact that we considered temperatures far 
beyond the range where TDGL is justified by microscopic theory. 

Nonuniformity leads to smaller persistent currents (in comparison with a uniform ring 
with the same average parameters), to smaller slopes of the current as a function of the fiux 
when passing between fiuxoid states, and to a larger temperature range for which the passage 
between fiuxoid states is non-hysteretic. It also leads to a narrower range of met ast ability. 
Finite width leads to smaller persistent currents (in comparison with ID rings) and to lesser 
sensitivity to nonuniformity. 

For several experiments in which the rings were intended to be uniform, better agreement 
with theory can be obtained by assuming that the rings actually had some unintentional 
nonuniformity; for those rings that were intentionally nonuniform, good agreement with 
theory was found using the reported values for the shape of the ring. 

Several years ago it was predicted that the passage between fiuxoid states may be medi- 
ated by a fiux-induced vortex^^ and this prediction has been experimentally confirmed.—*^ 
Our present study confirms the existence and position of this vortex. Moreover, it describes 
a situation in which the intermediate state between two fiuxoid states is actually a dynamic 
situation with sporadic migrations between both states, involving a place where the order 
parameter vanishes. 

The Kramer-Watts- Tobin (KWT) model was applied to two situations, mainly in order 
to estimate the expected discrepancy incurred when using TDGL. The first case was that of 
the persistent current in a fiuctuation range, i.e., a range for which the current would vanish 
in the absence of thermal fiuctuations. Remarkably, even for the large value of Tph in the 
case of aluminum, the currents obtained for KWT agree with those obtained for TDGL, and 
they both agree with the measured values. This agreement may be understood if we bear 
in mind that the persistent current is an equilibrium quantity and therefore depends only 
on the energy spectrum and not on the dynamics; since the energy functional for KWT is 
the same as for TDGL, they both lead to the same currents. In more general terms, we can 
say that even if some features of a model, such as TDGL, are not an accurate description of 
reality, the model can accurately describe situations in which these features are irrelevant. 

The second case in which we used KWT was the study of the lifetime of metastable states. 
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We find that the ring can remain in the metastable state for longer times than predicted by 
TDGL. Since according to Eqs. f|T2|) and f|T3|) KWT brings about a slower evolution of the 
absolute value of the order parameter and a faster evolution of its phase, we can conclude 
that the evolution rate of the absolute value is dominant in determining met ast ability. This 
can be understood, since the winding number of the order parameter has to change in order 
to escape from a fiuxoid state and this cannot happen as long as the absolute value of the 
order parameter is positive everywhere; on the other hand, when the order parameter does 
vanish at some point, the fast evolution of the phase could lead to multiple transitions, as 



discussed in Ref. 
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Appendix A: Numerical aspects of KWT 

In the case of TDGL, we may regard the evolution of the order parameter as isotropic in 
ip-space. On the other hand, Eqs. (fT2|) and (1131) show that according to KWT variations of 
ipk in the "radial" direction are inhibited by a factor h\^\, whereas in the "angular" direction 
they are enhanced by a factor For situations in which Klip^l'^ ^ 1, this factor can 

be important, and this forces us to divide steps in the angular direction into several steps, 
i.e., for every step lasting time r in the radial direction there are steps in the angular 
direction, each lasting time r/n^ {n^ ~ ^h)- 

In cases in which we evaluate an equilibrium property, it might be unnecessary to perform 
all the steps in the angular direction, since, due to the fast variation of Xk while \ipk \ remains 
frozen, we may assume that after a relatively small number of steps rip^oh ^ have 
already obtained the equilibrium value of the evaluated property for the given value of \ipk\, 
and we may go on and probe new values of \ipk\- In this case it is important to note that 
the statistical weight of this partial average is r and not Upj-ohT/n^. 
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FIG. 14: Size of the current as a function of temperature for (p = 0.25. The continuous hues 
were calculated using Eq. ()B7p and the values of I and Tc are marked next to each line; the other 

3^ . The mean free path obtained from normal conductivity 



parameters were taken from Ref. 
measurements is £ = 40 nm. The dashed line is the fit to the measurements given in Ref. 
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Equation f|T^ was implemented through the assignment ipkit+r) = ipk{t){'^+^\'^k\/\'^k\)- 
Some caution is required due to the presence of the factor \ipk\'^ in the denominator, which 
might accidentally be very small. We replaced it with \ipk\'^ + tiIj'^ /f; unless \ipk\ is exception- 
ally small, this replacement gives rise to a negligible O(r^) contribution, while it guarantees 
that ipk remains conveniently bounded in these exceptional cases. 



Appendix B: Persistent Currents in Normal Rings 

In normal metals current is due to single electrons, so that persistent currents are expected 
to have a 2$o flux period; the size and sign of the current are governed by coslkpL), which 
is random in practice. Recent experiment s22,i21 on single rings are consistent with theoretical 
predictions. 

On the other hand, several experiments on arrays of isolated rings of normal metals^"— 
exhibit periodicity, dependence on the number of rings, sign and orders of magnitude that 
would be appropriate for superconductors. Available theories^^"— indeed regard the average 
persistent currents as fluctuation superconductivity. 

Since one of the objectives of this article is to check the applicability of the Ginzburg- 
Landau model to situations that are not guaranteed by a microscopic treatment, we compare 
this model with one of the experiments. By using this model, we are considering thermal 
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fluctuations only, whereas quantum fluctuations are being ignored. According to Ref. |37|, this 
approach should give the dominant contribution to the current for temperatures > Tlp. The 
critical temperature will be regarded as a phenomenologic parameter, such that Tc > 
describes a superconductor, Tc = describes a perfect normal conductor, and the possibility 
of an "antisuperconductor" with < (which could be due to repulsive interaction among 
the electrons) can also be contemplated within this model. 

We consider a two-dimensional narrow uniform ring, write the Ginzburg-Landau energy 
in terms of the gauge-invariant order parameter, and neglect the quartic term. The energy 
becomes 



G 



h w 



2mD Jo ./n V 



D 



ds / dr 
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i^r + ivv'i- 
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with 1] = = 0.19kB{T — Tc)L'^/Mvf- ip obeys the boundary conditions iIj{s + 

L) = e'^'^^'^ip{s), where = $/$o, and dip/dr = at r = and r = D. 
Performing the transformation 
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The persistent current can be written as (/) = — (2ehw / mLD) (Im ds dr ilj*d'4>/ds). 
Using transformation ( 1B2I) we obtain 
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and introducing the statistical averages the current becomes 
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Summation over j gives 
(/) = - 
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with f{x) = xcothx; for the relevant experimental values, f{x) may be replaced with 1, 
which is the zero- width limit. 

The series in Eq. (IB6p is not absolutely convergent and we therefore need some physical 
criterion in order to sum it. Since n + (p is proportional to the velocity of Cooper pairs, 
and the sum should actually not contain the terms with velocities above the pair-breaking 
velocity, we introduce a cutoff for |n + 0|. This is done by replacing Eq. (]B6|) with 

/A = V (^ + 0)[l + /(2vrZ^V^ + (n + '/')V^)] .^7. 

^' nh ^ T] + (n + 4>y + e(n + ' ^ ^ 

n=— 00 ^ ' ^ ' 

where e is a small number. Taking e = 0.01 or e = 0.001, or adding e(n + </))^ rather than 
e(n + 0)"^, we obtain practically the same results. 

As a function of the flux, Eq. flB7p predicts a persistent current that vanishes when 
is an integer multiple of 0.5 and has maximum size when is close to an odd multiple of 
0.25. Figure [H] compares the size of the current at = 0.25 predicted by Eq. flB7p with 
the experimental values obtained for silver rings. The values of the parameters were taken 
from Ref. [S^. This experiment was not considered in Ref. |37|, because measurements were 



performed at a high frequency. In order to allow for the possibility that only a subset of 
the electrons on the Fermi surface are involved in superconductivity, we also considered 
values of the mean free path that are different from the average value deduced from normal 
conductivity. The best agreement with the experimental results is obtained for a rather 
large positive Tc, favoring the scenario advanced in Ref. |36 . 
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